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ABSTRACT 

The light curves of many Blazhko stars exhibit intervals in which successive pulsa¬ 
tion maxima alternate between two levels in a way that is characteristic of period¬ 
doubling. In addition, hydrocode models of these stars have clearly demonstrated 
period-doubling bifurcations. As a result, it is now generally accepted that these stars 
do indeed exhibit period-doubling. Here we present strong evidence that this assump¬ 
tion is incorrect. The alternating maxima likely result from the presence of one or 
more near-resonant modes which appear in the stellar spectra and are slightly but sig¬ 
nificantly offset from 3/2 times the fundamental frequency. We show that a previously 
proposed explanation for the presence of these peaks is inadequate. The phase-slip of 
the dominant near-resonant peak in RR Lyr is shown to be fully correlated with the 
parity of the observed alternations, providing further strong evidence that the process 
is nonresonant and cannot be characterized as period-doubling. The dominant near¬ 
resonant peak in V808 Cyg has side-peaks spaced at twice the Blazhko frequency. This 
apparent modulation indicates that the peak corresponds to a vibrational mode and 
also adds strong support to the beating-modes model of the Blazhko effect which can 
account for the doubled frequency. The modulation also demonstrates the “environ¬ 
ment” altering effect of large amplitude modes which is shown to be consistent with 
the amplitude equation formalism. 

Key words: instabilities - stars: oscillations (including pulsations) - stars: variables: 
RR Lyrae - stars: individual: RR Lyr - stars: individual: V808 Cyg 


1 INTRODUCTION 

It is well known that some nonlinear dynamical systems 
can e xhibit a period-doubling bifurcation (see, e.g. IStroeat j 
1 1994) . For parameter values below the bifurcation point, the 
system is oscillating periodically at frequency / and may also 
have spectral peaks at the harmonics of this frequency. In 
passing the bifurcation point, a small peak at //2 emerges 
and grows in size. The oscillations exhibited by the system 
alternate between two slightly different paths in the phase 
space, so that the period has suddenly become twice as long 
as it was originally. Typically this is observed in the suc¬ 
cessive maxima of the oscillation which alternate between 
two slightly different levels. Odd harmonics of //2 typically 
appear in the spectrum as well. We will refer to these as 
half-integer frequencies. Often additional period-doublings 
occur as the parameter is further advanced resulting in the 
appearance of //4, //8, etc. This cascade typically goes all 
the way to //oo over a Hnite change of the parameter and 


* E-mail: pbryant@ucsd.edu 


beyond this point the dynamics do not repeat and the sys¬ 
tem is called “chaotic”. 

iKolenberg et al.l ll2010l) hrst reported seeing alter¬ 
nating maxima in t he li ght curves of some Blazhko 
stars and lSzabo et al.l ll2010l) demonstrated period-doubling 
in a hydrocode model of RR Lyr. Since then it has 
been generally accepted that the observed alternation ef¬ 
fect was ind eed caused by a period-doubling bifurca¬ 
tion (see, e.g., iBuchler &: Kollathll201ll: iKollath ^ al.ll201ll: 


G_u^Eenberger_et_^ 20l4^ Mohiaxet^4 2012 |_|Kolenbe^ 


20l4 : lBenko et al.ll2ni4 : ISzabo et al. 12014 : iLe Borgne et al 
2014) . In this paper we show the serious problems with this 
claim, and present strong evidence to support an alternate 
explanation, namely that the alternation is caused by the 
presence of one or more excited modes whose frequencies are 
close to 3/2 times the fundamental frequency and (in most 
cases at least) are not in resonance with the fundamental 
mode. The largest of these peaks have frequencies that can 
be accurately determined and can be observed over the en¬ 
tire Kepler data set of about four years. Other studies have 
noticed that these peaks are off resonanc e, but a seemingly 
plausible explanation has been offered (ISzabo et al.l l201Cll . 
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Section 3.2) and as a result the idea persists that presence of 
these near-resonant peaks indicates that a period-doubling 
bifurcation has occurred. This belief is strengthened by 
the fact that actua l period-doubling h as been observed in 
hydrocode rn odels (ISzabo et al.l I2OI0I : iKollath et al.l I2OIII : 
ISmolec et al.ll2()llh . and it is assumed that therefore the 
same must be occurring in actual RR Lyrae stars. But as is 
shown below, these hydrocode results deviate quite substan¬ 
tially from the results obtained from actual stellar data and 
thus they do not prove that period-doubling is occurring in 
the real stars. _ _ 

The argument by ISzabo et al.l ll201(t) is that the half¬ 
integer frequencies appear to be off-resonance because they 
are effectively modulated, periodically and/or randomly, 
and that this produces a cluster of side-peaks in the spec¬ 
trum. It does not explain however, how the central peak can 
be shifted to one side of the exact half-integer frequency. 
Even for a completely random modulation the result should 
be a cluster of peaks that is centered on the base frequency. 
This is in fact what happens in their simulation of the effect 
as seen in their Figure 10, where the cluster of peaks in part 
(b) does line up with the main peak in part (a), and if the 
image is magnified it can be seen that the main peak in part 
(a) is in precise alignment with one of the peaks in part (b), 
and therefore this peak is exactly on-resonance. (Of the two 
tallest peaks in the cluster, it is the one on the left.) 

We will show that the dominant peak near (3/2)/o in 
RR Lyr is substantially off-resonance and persists over the 
entire Kepler data set and is therefore not an artifact of 
some random or chaotic process. We will also show that 
the observed half-integer oscillation is not phase locked to 
the fundamental but is constantly slipping in phase at the 
rate predicted by the offset of the peak location from ex¬ 
act resonance. We also show that a similar off-resonance 
peak exists in the spectrum for V808 Cyg, another star 
whic h is believed to e xhibit a strong period-doubling ef¬ 
fect dSzabo et al.ll20i^ . This peak has side-peaks consis¬ 
tent with modulation at double the Blazhko frequency. One 
explanation for his strange doubling effect (described in de¬ 
tail below) is that this peak does in fact correspond to a 
nonresonant mode and that the side peaks result from in¬ 
teractions between this mode and the modes that gener¬ 
ate the Blazhko effec t as described by the “beating-modes 
model” (lBrvantll2015l )Fl It thus simultaneously presents ev¬ 
idence against period-doubling and evidence in favor of the 
beating-modes Blazhko model (since it would appear that 
other models cannot provide an explanation for the doubled 
frequency). The analysis introduces the concept of the “en¬ 
vironment” within which the modulated mode resides and 
which is generated by the larger amplitude Blazhko modes. 
We show that this c oncept is compatible with the amplitude 
equation formalism llBuchler Sz Goupillll984f l. 

We note that the presence of a near-resonant mode in 
the beating-modes model is very similar to the currently 
proposed explanation for the maxima alternation, but aside 
from this similarity the two explanations are not linked, i.e. 
the validity of the current work does not depend on that 


^ Not to be confused with iBrvand l |201^, an earlier work by the 
author on the Blazhko effect, which iBrvantI (I 2 OI 5 I I reinterprets 
and supersedes. 



Figure 1. Spectra for RR Lyr in the vicinity of (3/2)/o computed 
from Kepler project KIC 7198959, long cadence corrected flux 
data using Period04 software. Average (or zero point) flux for 
this data is about 1.2222 X 10^. Data was pre-whitened, removing 
the fundamental and its harmonics through tenth. The upper 
curve is computed from quarters Q5, Q6 and Q7; it is offset by 12 
units vertically for clarity. The second curve uses Q7, Q8 and Q9 
with offset 9. The third curve uses Qll, Q12 and Q13 with offset 
6. The fourth curve uses Q14, Q15 and Q16 with offset 3. The 
fifth curve uses all available data: Q1 through Q17 with no offset. 
Note that Q1 and Q17 are partial quarters and Q3, Q4 and QIO 
are not available. The dashed vertical line marks the location of 
(3/2)/o. Note that none of the five curves have a peak at this 
location, while they all have one at about 2.6645 which is clearly 
nonresonant. 

of the Blazhko explanation or vice-versa. Note that in the 
beating-modes model there can be nonlinear interactions 
between these two modes as long as these interactions do 
not cause phase-locking. Some nonlinear interactions will 
depend on the relative phase of the two modes, but since 
this is constantly changing, it will typically have a modu¬ 
latory effect, going through one modulation cycle with the 
periodicity of the relative phase. But this modulation is a 
secondary effect, the primary cause of the Blazhko effect is 
the phase-slip itself. iBrvantI (|2015l ) reproduced with good 
accuracy the Blazhko effect of RR Lyr without including 
any actual modulation. One might assume then that such 
modulation effects are relatively small for RR Lyr, but this 
is not necessarily the case for other Blazhko stars. 


2 RESULTS AND DISCUSSION 

In Figure[T]we show spectra for RR Lyr. The first four curves 


© 2015 RAS, MNRAS 000, [Hg] 




































Is there period-doubling in Blazhko stars? 3 


Table 1. Correlation between the phase slip of the near resonant 
mode and the phase of the maxima alternation pattern for RR 
Lyr. A single point is chosen near the middle of each low-high 
alternation sequence. It is required to be the point of maximum 
flux in a “high” cycle. In some cases the sequence may include a 
few bad or non-alternating points provided the alternation con¬ 
tinues with the proper phase. The first point in the table has 
truncated Barycentric Julian Date 54992.569. Columns: time is 
the time in days relative to the first point; len is the approximate 
length in cycles of the alternation sequence; eye is the cycle num¬ 
ber, calculated by multiplying time by fo (1.76429291); pari is 
the parity of eye after rounding to the nearest integer; slip is the 
phase slip in half cycles, calculated by multiplying time by twice 
the frequency offset of the near resonant mode (2 x 0.0181) and 
adding a correction of -0.3 (chosen to make the best overall fit 
to the entire set); par2 is the parity of slip after rounding to the 
nearest integer; cor is the correlation of pari and par2, indicated 
as “yes” if they match and “no” if they do not. 


time 

len 

eye 

pari 

slip 

par2 

cor 

0 

12 

0 

even 

- 0.3 

even 

yes 

36.843 

24 

65.001 

odd 

1.034 

odd 

yes 

65.715 

16 

115.941 

even 

2.079 

even 

yes 

285.658 

12 

503.985 

even 

10.041 

even 

yes 

306.644 

54 

541.011 

odd 

10.801 

odd 

yes 

338.931 

24 

597.973 

even 

11.969 

even 

yes 

392.244 

32 

692.033 

even 

13.899 

even 

yes 

424.55 

30 

749.03 

odd 

15.069 

odd 

yes 

456.814 

32 

805.954 

even 

16.237 

even 

yes 

540.159 

18 

953 

odd 

19.254 

odd 

yes 

587.789 

11 

1037.032 

odd 

20.978 

odd 

yes 

618.95 

35 

1092.009 

even 

22.106 

even 

yes 

653.504 

26 

1152.972 

odd 

23.357 

odd 

yes 

669.974 

17 

1182.03 

even 

23.953 

even 

yes 

701.136 

20 

1237.01 

odd 

25.081 

odd 

yes 

731.175 

43 

1290.007 

even 

26.169 

even 

yes 

926.722 

14 

1635.009 

odd 

33.247 

odd 

yes 

1036.675 

20 

1828.998 

odd 

37.228 

odd 

yes 


from the top down were calculated for groups of three con¬ 
secutive quarters of Kepler data as indicated in the caption. 
The bottom curve was calculated using all available Kepler 
data, quarters 1 through 17, and thus has much higher fre¬ 
quency resolution. Note the consistent pattern in this set of 
curves with a dominant peak at 2.6645, henceforth called 
/a, while the expected location for resonance would be at 
2.6464 (marked with a dashed vertical line in the figure). 
We propose that /a is the base frequency of a mode which 
we will refer to as Mode A. The frequency difference being 
0.0181 means that this oscillation will slip in phase by one 
half cycle about every 27.6 days. This phase slip can be ob¬ 
served in the alternation pattern, since a half cycle slip will 
cause a switch fr om the even cycles bei ng high to the odd 
cycles being high. lLe Borgne et al. I ||2014| also found this fre¬ 
quency (0.018) in their spectral analysis of the alternation 
effect (see their Figure 3). We verify this phase slip effect in 
Table [T] which shows a full correlation in the data for RR 
Lyr between the high cycles being odd or even and the num¬ 
ber of half cycle slips being odd or even. Most, if not all, of 
the easily visible alternation sequences of length longer than 
10 are included in the table. Note that this continuous phase 
slip means that there is no phase-locking between the funda¬ 
mental and Mode A and therefore there is no resonance and 


no period-doubling. The fact that this effect can be followed 
through the entire 17 quarters of Kepler data indicates that 
this is tied to something with a very stable frequency (i.e. 
a non-resonant mode). It can’t be attributed to a random 
phase drift and it seems quite incompatible with period dou¬ 
bling. 

Note that the alternation amplitude goes to zero at the 
transition point between where the high cycles are even and 
where they are odd. This happens again when they transi¬ 
tion back to even, and so on. The amplitude of Mode A has 
not gone to zero at these points, but its relative phase is 
such that it is not visible as a maxima alternation. 

An interesting thing to note in Figure[T]is that the peak 
/a is diminishing in height over time. The amplitude of the 
Blazhko effect is also decreasing simultaneously. This might 
be taken as evidence of a possible connec tion between t hem, 
but it could also be just a coincidence. iBrvand ll2015ll has 
shown that the Blazhko effect may also be an artifact of 
a near-resonant nonradial mode. So one simple explanation 
for the correlation could be that the excitation levels of all 
of these nonradial modes are decreasing for some reason. For 
example, if the magnetic field of the star was slowly increas¬ 
ing in time, this might cause a slow increase in turbulence 
resulting in increased damping and/or decreased excitation 
of the modes and this would affect both the Blazhko ef¬ 
fect and the maxima alternation. A similar magn e tic me ch- 
anism was recently proposed by IStotherj ll2006l . I 2 OI 0 II as 
part of an explanati on of the Blazhko e ffect. One critique 
of this explanation llMolnar et al.l l2012l l claimed that the 
mechanism was too slow to explain the Blazhko effect, but 
this would not disqualify it from explaining these long term 
fluctuations in the Blazhko effect. But whatever the reason, 
long term variations in the Blazhko effect in RR Lyr have 
been known for some t ime and remain unexplained; see e.g. 
iLe Borgne et akl (l2014h . Period doubling does not provide a 
definitive explanation for this mystery either. The observed 
correlation between the diminishing of /a and the dimin¬ 
ishing of the Blazhko effect could be real or it could be a 
coincidence, there is insufficient information to decide based 
on data coming from just one star. However, if it was real, 
it could explain why the maxima alternations have thus far 
only been found in Blazhko stars. 

The spectrum also shows additional peaks at frequen¬ 
cies that are mixing products of this peak with the fun¬ 
damental and its harmonics, i.e. frequencies of the form 
/a -b nfo where n is an integer and fo is the fundamental fre¬ 
quency. These peaks are all near, but not equal to, the half¬ 
integer frequencies that are expected for period-doubling. 
The amplitudes of these peaks are determined, based on 
all available data, and the results compared in Figure [2] 
with the corresponding spectrum for period-doubling in a 
hydrocode model o f RR Lyr (taken from M8 in Figure 11 in 
ISmolec et al.ll2()lll L Note that their spectrum has the domi¬ 
nant peak at (l/2)/o rather than at (3/2)/o. The ratio of the 
peak heights (the 3/2 peak height over the 1/2 peak height) 
changes rather drastically between the two cases. This lack 
of agreement in the spectrum can be taken as fairly strong 
evidence that the model is not presenting a correct represen¬ 
tation of the dynamics of the actual star. Having the max¬ 
imum half-integer peak near 3/2 in the observational spec¬ 
trum is, of course, consistent with the source for this peak 
being an excited mode with that frequency. The other half- 
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Figure 2. Comparison of the observed half-integer frequency 
peaks with those of a hydrocode model exhibiting period¬ 
doubling. The solid curve with diamond data points is calculated 
from Kepler data for RR Lyr using all available quarters, Q1 
through Q17. The vertical scale for this curve is on the left. These 
peaks are off-resonance as discussed in the text. The dashed curve 
with inverted tr iangle data po i nts is replotted from the M8 curve 
in Figure 11 in ISmolec et al.l ll201lh . The vertical scale for this 
curve is on the right. These peaks are exactly on-resonance. Note 
the great disparity between the two curves, the first having a 
strong maxima at 3/2 while the second has a strong maxima at 
1 / 2 . 

integer peaks are then nonlinear mixing peaks, which would 
typically be expected to be smaller in amplitude than the 
one near 3/2 in agreement with observations. Other obser- 
vatio nal results also ha ve the maximu m at 3/2 rather than 
1/2 (ISzabo et al.ll201ol . Figure 8). The lSmolec et ahl (1201 il l 
results also show an exact correlation between the Blazhko 
phase of their model and the appearance and disappearance 
of the period-doubling, something that is entirely missing in 
the actual stellar data. This is of course only relevant if we 
assume that the Blazhko effect involves actual modu lation 
as opposed to a beatin g-modes process (lBrvantl[2015l l. 

As was noted by ISzabo et al.l (l201ol l. when observed 
on a short timescale the amplitudes of the half-integer fre¬ 
quency peaks appear to be modulated in an apparently ran¬ 
dom fashion (see their Figure 3, bottom panel). Here are 
three possible reasons (all of which could be acting simul¬ 
taneously): 1) Turbulence noise is interacting with Mode 
A and causing these fluctuations. Since this mode is only 
marginally unstable, it wil l be very sens i tive and can ef¬ 
fectively amplify this noise. I^uchler et al.1 (Il993l 'l previously 
studied such stochastic excitation of stellar pulsations. 2) 
Some of the nearby peaks in the spectrum may correspond 
to separate excited modes and add additional near-resonant 
frequencies into the system. Several of these added together 
could generate something that looks like a random modu¬ 
lation. 3) Noise and measurement error of the Kepler elec¬ 
tronics could be affecting the results. 

This modulation or apparent modulation of Mode A 
also means that occasionally the observed maxima alterna¬ 
tions are considerably larger than might be expected from 
its peak height in the spectrum. But since the time series 
can always be reconstructed from the spectral information, 
these occasional large alternations must be reflected in the 
combination of the main peak with various associated side 
peaks and harmonics when these are all in phase. 

An interesting question is why /a is so close in fre¬ 
quency to (3/2)/o. The usual mechanism for transferring 



Figure 3. Spectra for V808 Cyg in the vicinity of (3/2)/o com¬ 
puted from Kepler project KIC 4484128, long cadence corrected 
flux data using Period04 software. Average (or zero point) flux for 
this data is about 8.8363 X 10®. Data was pre-whitened, remov¬ 
ing the fundamental and its harmonics through tenth. The upper 
curve is computed from quarters Q2, Q3 and Q4; it is offset by 
15 units vertically for clarity. The second curve uses Q7, Q8 and 
Q9 with offset 10. The third curve uses Qll, Q12 and Q13 with 
offset 5. The fourth curve uses all available data: Q1 through Q17 
with no offset. Note that Q1 and Q17 are partial quarters and Q6, 
QIO and Q14 are not available. The dashed vertical line marks 
the location of (3/2 )/q. Note that none of the five curves have 
a peak at this location, while they all have one at about 2.6983 
which is clearly nonresonant. See text for discussion of the split 
peak in the bottom curve. 

energy from the fundamental mode to another mode only 
works when there is phase-locking between the modes and 
that is apparently not happening in this case. Bnt perhaps 
there is some other, as yet unknown, mechanism at work. It 
is also possible that the proximity is merely an accident. In 
Figure [T] note the two smaller peaks to the left of the main 
peak. There are more peaks further to the left (off scale) of 
a similar size and more densely packed together. This sug¬ 
gests that there is a cluster of peaks that just happens to 
slightly overlap the resonance point. Other stars, most no¬ 
tably V445 Lyr, appear to hav e lots of excited modes tha t 
are far from resonance; see e.g. iGuggenberger et al.1 ll2012l) . 
If a mode can be excited in arbitrary locations then near¬ 
resonant locations should also be possible. 

In Figure [3] we show spectra for V808 Cyg. The first 
three curves from the top down were calculated for groups 
of three consecutive quarters of Kepler data as indicated 
in the caption. The bottom curve was calculated using all 
available Kepler data, and thus has much higher frequency 
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resolution. Note the consistent pattern in this set of curves, 
with a dominant peak at 2.6983, which we will again call 
/a, while the expected location for resonance would be at 
2.7379 (marked with a dashed vertical line in the figure). 
This time the dominant peak is found on the opposite side 
of the resonance compared to Figure [T] An interesting fea¬ 
ture of the bottom curve is that the /a has split into two 
peaks with frequencies of 2.6977 and 2.6999. One explana¬ 
tion would be that the noisy turbulent environment of the 
star is causing the amplitude to fluctuate in a way that hap¬ 
pens to look like the sum of these two frequencies. Another 
explanation is that there really are two closely spaced modes 
that generate these peaks. In favor of the second explanation 
is the fact that the difference between the peak frequencies is 
0.0022 and from this one can determine that the two modes 
slipped in phase by about 3.2 cycles over the time span of 
the full data set. This seems like an excessive number of 
consistent cycles to be accidentally generated by a random 
process. A second interesting feature is the fact that the pair 
of peaks appears to have pairs of side peaks. A total of 8 side 
peaks (four pairs) are visible (although one is barely visible 
above the noise). The fact that it appears modulated is a 
strong indication that this peak (or pair of peaks) does cor¬ 
respond to an actual vibrational mode (or pair of modes) 
as opposed to the explanation that it is some kind of noisy 
period-doubling byproduct. 

Perhaps even more interesting is that the apparent mod¬ 
ulation frequency, based on the side-peak spacing, is twice 
the Blazhko frequency. By our measurements from the en¬ 
tire data set, the Blazhko frequency is about 0.01085, so 
twice this is 0.02170. We measure the side-peak spacing 
to be about 0.02165 so this is a very accurate match. We 
have an explanation for this doubled frequency i n terms 
of th e beating-modes model of the Blazhko effect dBrvantl 
I 2 OI 5 II . That model proposes that the apparent modulation 
of the fundamental is actually the beat frequency produced 
by simply adding the fundamental mode oscillation to that 
of a near-resonant (presumably nonradial) mode, whose fre¬ 
quency differs from the fundamental by the Blazhko fre¬ 
quency. As we explain below, the modulation of Mode A 
observed in the current problem may be due to the type of 
nonlinear interaction that is not phase-dependent. 

Let us assume, to be specific, that the nonradial mode 
is a dipolar mode with axial s ymmetry, i.e. a mo de with in¬ 
dices I — 1 and m = 0 (see, e.g. lUnno et al.lll989l . Section 4). 
The dipolar character will cause the radial component of the 
oscillations in the northern hemisphere to be be 180 degrees 
out of phase with those in the southern hemisphere. If, at 
one point in time, this oscillation has a particular phase re¬ 
lationship with the fundamental mode in the northern hemi¬ 
sphere, then one half of a Blazhko period later it have the 
identical phase relationship with the fundamental mode in 
the southern hemisphere. So for Mode A, these two points 
in time present a spatially inverted but otherwise equivalent 
“environment” for that mode to reside in. Thus if there is 
any variation in this environment then it must repeat every 
half Blazhko cycle, i.e. it must oscillate at twice the Blazhko 
frequency. 

This type of oscillatory environment can have a modu¬ 
latory effect on the frequency of Mode A. This can be seen 
through the cross coupling “T” terms tha t appear in the am- 
plitude equations, see, e.g. Equation fll lBuchler &: KollathI 


(2011), Equations (la) and (lb) in |Moskalik_&_^uchlet 


199C) and Equations (70) and (71) in iBuchler fc Goupil 


198411 . In the equation for da/dt, this term has the form 


T\bfa where T is a complex coefficient and a and b are 
the complex amplitudes of Modes A and B. The imaginary 
component of T results in the frequency of Mode A being 
shifted in proportion to the square of the amplitude of Mode 
B, while the real component of T results in a similar shift in 
the growth rate of Mode A. Note that this term, unlike some 
coupling terms, is a non-resonant term. By this we mean 
that it does not depend on the relative phase of Mode A and 
will not cause phase-locking between the modes. This is in 
contras t, for example, to the “C” terms appearing in Equa¬ 
tion fll lBuchler fc KollathI (l 201 lll which can result in phase¬ 
locking between modes which might otherwise be slightly off 
resonance. So the presence of the T term in the equation for 
da/dt can be thought of as an effect of the “environment” 
generated by Mode B on the dynamics of Mode A. 

Extending this analysis to the current problem with 
three modes, we will let Mode B be the fundamental mode, 
and let Mode C be the mode that is near-resonant with the 
fundamental (according to the beating-modes model), while 
Mode A retains its definition from earlier in the paper. To 
represent a combined effect of modes B and C in generating 
an environment for mode A, we may add {U\b*c + U 2 C*b)a 
to the equation for da/dt, where U\ and U 2 are complex 
coefficients. Unlike the T terms, these generate an environ¬ 
ment that is slowly varying and will therefore slowly change 
the frequency and/or amplitude of Mode A. Since the fre¬ 
quency difference between Modes B and C is the Blazhko 
frequency, the result will be to modulate Mode A at that 
frequency. If we assume that mode C is a dipolar mode, then 
one might expect, as a result of the environmental symmetry 
arguments given above, that the coefficients U\ and U 2 for 
this case would be exactly zero. However these arguments 
would not apply for axisymmetric modes with even sym¬ 
metry about the equator, i.e. those with even degree I and 
order m = 0. So for those cases these terms could produce 
modulation of mode A at the Blazhko frequency. For axisym¬ 
metric modes of odd degree I (including the dipolar modes) 
the modulation frequency should be twice the Blazhko fre¬ 
quency and the lowest order environmental terms that could 
produce this result would appear to be of the form (Vib*^c^ 
+V 2 C*^ 6 ^)a. Note that as the degree I increases, this distinc¬ 
tion between odd and even diminishes in importance and 
so the U coefficients for even I should tend toward zero. 
Also note that modes with nonzero order m cannot gener¬ 
ate this type of oscillatory environment in conjunction with 
the fundamental unless they are excited in pairs with whose 
m values add to zero. In this case there would again be a 
symmetry that would only allow modulation at twice the 
Blazhko frequency. 

An interesting question is whether the symmetry prop¬ 
erties of the dynamics are retained when the amplitude is 
large and therefore strongly nonlinear. This may be indi¬ 
cated in the current case by the highly non-sinusoidal nature 
of the fundamental oscillation. Studies of simpler nonlin¬ 
ear systems have indicated that symmetrical quasiperiodic 
oscillations may be qui te common; see e .g. Figure (2b) in 
iBrvant Sz Jeffriea (1 19871 1. As discussed in iBrvant Sz JeffriesI 
( 19841 . 1987 1 another possibility is complementary asym¬ 
metric attractor pairs, though this would seem unlikely ex- 
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cept when phase locki ng occurs; see e.g. Figure (3b) in 
iBrvant &: Jeffriej dlQSifl . 
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